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Abstract 

We study the eigenvalue distribution of a random matrix, at a transition where a 
new connected component of the eigenvalue density support appears away from other 
connected components. Unlike previously studied critical points, which correspond to 
rational singularities p(x) ~ x p ^ q classified by conformal minimal models and integrable 



hierarchies, this transition shows logarithmic and non-analytical behaviours. There is 
no critical exponent, instead, the power of N changes in a saw teeth behaviour. 

1 Introduction 

Random matrix models (201 EH! have been studied in relationship with many areas of 
physics and mathematics. The reason of their success for most of their applications is 
their "universality" property, i.e. the fact that the eigenvalues statistical distribution 
of a large random matrix depends only on the symmetries of the matrix ensemble, and 
not on the detailed Boltzmann weight (characterized by a potential). Although this 
universality property has been much studied for generic potentials, some universality 
should also hold for critical potentials. Different kinds of critical potentials have been 
studied, and their universality classes have been found to be in correspondence with 
non-linear integrable hierarchies (KdV, MKdV, KP,...) [3 13 HQ El E2| , and with 
the (p, q) rational minimal models of conformal field theory [§.• . They correspond to 
rational singularituies of the equilibrium density: 

p{x) ~ (x - a) p/q . (1.1) 



1 E-mail: eynard@spht.saclay.cea.fr 
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In the hermitian 1-matrix model, we have only q = 2 and p arbitrary, thus we get 
only half integer singularities (hyperelliptical curves), which are in relationship with 
the KdV hierarchies, whereas a 2-matrix model allows to have any rational singularity 
(p, q) 8j. The specific heat near such a rational singularity obeys a Gelfand-Dikii-type 
equation (Painleve I equation for (p,q) = (3,2)). A well known case is the edge of 
the spectrum where (p,q) = (1,2), which gives Tracy- Widom law [2H1, and which is 
governed by the Painleve II equation. Another well known case is the merging of two 
cuts (Bleher and Its 0E]), where (p, q) = (2, 1), which is also governed by a Painleve 
II equation as well (indeed, (p, q) and (q,p) are known to be dual to each other [Ej). 

Here, we shall study a kind of critical point which has been mostly disregarded 
(because usual methods don't apply to it): "the birth of a cut critical point" 2 . 

Such a critical point, is characterized by the fact that when a parameter of the 
model (let us call it temperature) is varied, a new connected component appears in the 
support of the large N average eigenvalue density. When the temperature T is just 
above critical temperature T c , the number of eigenvalues in the newborn connected 
component is small (see figHJ), and thus, many usual large n methods don't work in 
that case. 

Our goal is to study the eigenvalues statictics in the vicinity of the critical point, 
and find its universality class. 

In this purpose, we shall start from the partition function, and treat the eigenvalues 
in the other cuts with mean field approximations, and reduce the problem to an effective 
partition function for eigenvalues in the newborn cut only, in a method similar to j3j. 

We find that, unlike rational critical points, the birth of a cut critical point does 
not corrsepond to power law behaviors or transcandental differential equations, but it 
exhibits logarithmic behaviors, and discontinuous functions. 

The matrix model is associated to a family of orthogonal polynomials, whose zeroes 
lie inside the connected components of the density 1221 UH] • Our goal is also to study 
the asymptotic behavior of the orthogonal polynomials in the vicinity of the newborn 
cut. 

Outline of the article: 

- In section 2 we introduce definitions and notations for orthogonal polynomials and 
associated quantities. 

- In section 3 we recall classical results of random matrix theory: the semiclassical 
behaviors of free energy, density, correlation functions, orthogonal polynomials, valid 
away from critical points. 

2 Name suggested by P. Bleher who initiated this work. 
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Figure 1: The "birth of a cut" critical potential is such that one of the potential wells 
of the effective potential is just at the Fermi level. At a z/ th order critical point, the 
effective potential behaves like V e g(e) + ^V^ u \e) (x — e) 2u + . . .. 



- In section 4 we study the analytical continuation of the previous semiclassical ap- 
proximations, near the "birth of a cut" critical point (divergencies at critical point). 

- In section 5 we compute the partition function with mean field theory for eigenvalues 
not in the newborn cut, and derive an effective partition function for the newborncut 
eigenvalues. 

- In section 6 we use the results of section 5 to deduce the asymptotic behaviors of 
correlation functions and orthogonal polynomials in the vicinity of the critical point. 

- Section 7 is the conclusion. 

- In appendix, we recall Stirling's formula and its consequences, and we recall some 
elliptical function basics. 



Given an integer n (we will later consider the limit n — > oo), a real polynomial V(x) 
called the potential: 



2 Setting 
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Figure 2: The "birth of a cut" density of eigenvalues is such that one of the connected 
components of the support contains a very small number of eigenvalues (<C n). 



of even degree and positive leading coefficient (ga+i > 0), and a temperature T > 0, 
we define the partition function: 

„ n 

Z n (T, V):=- d Xl ... dx n (A(x t )) 2 IT (2.3) 
n. J R n . =1 

(where A(xj) = n^^i — x «) ^ s the Vandermonde determinant) and the free energy 
F n (T,V): 

e -^F n (T,V) ;= Z w (r,VQ (2 4) 

where iJ„ is a combinatorial normalization: 

n_1 3/4n 2 

H n := (27r)"/ 2 n"" 2 / 2 e 3 / 4 " 2 Y\ k\ = T^ 2 vr" 2 / 2 n^ 2 (2.5) 

and C/„ is the volume of the group U(n)/U(l) n . 
Then we define the resolvent: 



W n (x,T,V) :=—I—- [ dx 1 ...dx n —t— A 2 ( Xi ) f[e ~* 
n\Z n {T,V)J x-xi "t 



V(xi 



and its first moment: 



2m 

Notice that more generally for k > 0: 

dF n (T,V) 1 /■ , 



dg k 2m 



(2.6) 



r„(T, V) := n] - J y> y) y dn . . . da; n a* A 2 (a;,) f[e-^) 

i- y xW n (x, T, V)dx = V) (2-7) 



j x k W n (x,T,V)dx (2.8) 
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Then, given a temperature T c and an integer N, we define: 

IT "+1 T 
n+lK-l-c^-, v 

Zn(T e ±V) 



— v~FF~H TF\ > ^- y J 



and 



7 " :=l fc = V — — • (2 - 10) 

and 

«.==s(WT.=±i,v)-r.(T4,v))=-^. ( ,n) 

Notice that Z n (T, V), h n and 7„ are strictly positive for all n, in particular they don't 
vanish. 

We also introduce the functions [221 1201 UD] : 

Z n (l cW , V) y/h n 
which form an orthogonal family of monic polynomials (deg7r n = n): 

J 7r n (07r m (Oe _ ^ y(€) ^ = M™, , (2.13) 
and which satisfy the 3-terms recursion relation: 

£7T n (0 = 7T n+1 (0 + /W0 + 7>»-l(0 ■ (2-14) 

And we introduce the functions: 



, «« = 5#e#n« (, 15) 

which are the Hilbert transforms of the 7r n (x): 

= / -^^OOe-W (2.16) 
and which satisfy the same 3-terms recursion relation, with an initial term: 

£7Tn(0 = 7T n+1 (0 +^„7T„(0 +7n^n-l(0 + <*n,0 V (2-17) 

We also define the kernel: 

n— 1 -. 

y) ■= E ^ ^(OTi(y) e -&< v «> +v <«» , (2.18) 
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and it is well known that all correlation functions can be expressed in term of that 
kernel (Dyson's theorem JEJI^OJ): 

Pn{x) = -K n (x,x) 

n 

Pn(x, y) = -J (K n (x, x)K n (y, y) - K n (x, y)K n (y, x)) 
n 

(2.19) 

and that we have the Christoffel-Darboux theorem: 

K ( x y)= 1 MxW-ijy) - n n -i{x)ir n {y) e -£r(v(x)+v( v )) (2 20) 

Our goal is to study 7 n , j3 n , n n , n n , in the vicinity of N — > oo, and |n — <C AT, 
and with T c chosen such that we are at a special critical point described below. 
For the moment, let us study the large iV-limits away from the critical point. 

3 Classical limits 

It is well known that in the semiclassical limit iV — > oo, and \n — N\ <C N, if T ^ T c , 
the free energy has a large n limit [TTJ [THl EJ CHI UH] : 

F n (T, V) — > F(T, V) + 0(l/n 2 ) (3.21) 

and so has the resolvent W n (x, T, V): 

W n (x,T,V) -^W(x,T,V) + 0(l/n) (3.22) 

3.1 The large n resolvent 

It is well known that, if the potential is such that V'(x) is a rational fraction, with its 
poles outside the cuts (that assumption will become clear below), the large n resolvent 
W(x,T, V) can be written as the solution of an hyperelliptical equation [T3| IE]: 

W(x, T,V) = ± (V'(x) - M(x, T, V) a/ a(x, T, Vfj (3.23) 
where a is a monic even degree (2s > 2) polynomial with distinct simple zeroes only: 

s 

cr(x) = Y\( X - a i)( x - h) , . . . < di < bi < a i+1 < . . . (3.24) 

i=l 

whose zeroes are called the endpoints, and Ui=i[ a «> &i] i s called the support, and M is 
a rational function with the same poles as V . If one assumes that s and a are known, 
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M(x,T,V) is determined by the condition that W(x,T,V) is finite (in the physical 
sheet) when x —>■ oo and when x approaches the poles of V'(x). 

The large n limit of the density of eigenvalues is then: 



1 s 
p(x,T,V) = —M(x,T,VW-a(x,T, V) , xe\J[a h h] (3.25) 

8=1 

We also define the effective potential: 



V eS (x,T,V) :=V(x) -2Tlnx-2 / ( W(x, T, V) ) dx (3.26) 



x 



Notice that its derivative is V'{x) - 2W(x, T, V) = M(x, T, V)y/a(x,T,V). 

So far, we have not explained how to determine s and the polynomial a . If one 
assumes that s is known, a(x, T, V) is determined by the conditions: 



W(x,T,V) ~ - + 0(l/x 2 ) 

x-+co x 

if s > 1, Vi = 1, . . . , s - 1, V eS (bi) = V e $(a i+1 ) 



(3.27) 



The large n free energy is then given by: 



F(T, V) = -^j> W(x, T, V)V(x)dx + \TV eS {b s ) (3.28) 

where the integration contour is a counter clockwise circle around oo. 

The number of endpoints s = s(T, V), (we have 1 < s < d) is determined by the 
condition that the free energy is minimum (one can prove that s(T, V) < (d + l)/2) 

m 

3.2 Derivatives with respect to T 

Let us introduce: 

n{x ^ vy ^jn^,9^m (3 , 9) 

CJ ^o-(x, T, V) 

where Qn{x, T, V) is a monic polynomial of degree s — 1, determined by the conditions: 

if s > 1, Vi = l,...,s-1, ^t£^Hdx = (3.30) 

A, y/(r(x,T,V) 

In algebraic geometry, is called "normalized abelian differential of the third kind" 

naim. 

We introduce the multivalued function A(x,T, V): 

A(x, T, V) := exp Qf fi(x', T, V) rfx'^ (3.31) 



and 

7 (T,V):= lim * (3.32) 
ic-*oo A(x, i , l/J 



Then we have the following derivatives: 

V cS (b s ) (3.33) 



dF 
df 



d 2 F 
df 1 

dV eS (x) 
dT 



-2 In 7 (3.34) 
-21n( 7 A(x)) (3.35) 



Notice that: 



- = -|xO(x)dx (3.36) 
V eS (b s ) = ^- iW-2Tm 7 (3.37) 



3.3 Poles of the potential 

Assume that V'(x) has a simple pole at x — £, with residue r (it may have other poles 
too), then we define the function: 

C, T, V) := = ^= ^ — Q H (x, j (3.38) 

where Qh(x,0 is a monic polynomial in x, of degree s— 1, determined by the conditions: 

/■^Q ff (z,£,T,V) + ^ 
if s > 1, Vi = l,...,s-1, / — dx = 3.39 

We also define its (multivalued) primitive: 

hiE(x,g):= H{x',£)dx' (3.40) 

J oo 

Notice that it is finite near the endpoints and near x — £. In algebraic geometry, 
(x — £) / E(x, £) is related to the "prime form" [T%| fT7] . 
Then we have: 

()1,li, r) = ln(x-0 -21n£(x,0 (3.41) 



dr 
c 

~dT 2in 



dT 1 

= _ 4> x H(x, dx (3.42) 



^ = l(V(x)-V«(x))\ x= t (3.43) 
ln( 7 A(0) = ln(C - &,) - 21n£(6 s ,0 (3.44) 



drdT 

If V(x) has simple poles at x — £1 with residue ri and at x = £2 with residue r 2 , 
we have: 

d 2 F 

= ln^(6,6) (3-45) 



<9ri<9r 2 

and thus it is clear that lnE" has some symmetry properties: In E(x,y) = In E(y, x). 
3.4 1-cut case 

If W(x, T, V) has one cut [a(T, V), 6(T, V)] with a < 6, we use the Joukowski's param- 
eterization: 

b + a b — a 

x= — 1 — cosh<p (3.46) 

i.e. 

\/<t{x) = — ^ sinh (3.47) 

We have: 

Q(x) = - 1 = = 9 4 , A(s)=e«*> , 7=^ (3.48) 
V ; y/(x-a)(x-b) dx v ; 4 v ' 



H( x n - d ^ - e(x n = i- e -w x ) + *«» = x ~^ 

cfe e ^)+^)-l ' 1 ,U A(x)-A(0 

(3.49) 

Then, it is well known [] that we have the large n, N asymptotics (in the regime 
n/N =finite): 

b(T^) - a(T c f ) 6(r c $) + a(T c f ) 

7n ~ ^ , Pn ~ " {S.bL) 

3.5 2-cut case 

If W(x,T,V) has two cuts [a(T, V), b(T, V)] U [c(T,V),d(T,V)\ with a < 6 < c < d, 
Let m be their biratio: 

m = ^1 (3-52) 
(c — a) (a — 0) 
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We parameterize: 

d — a 



x(u) =d- - b _ a — (3.53) 



1 + ^fsn 2 (w,m) 



where sn is tha elliptical sine function (see appendix I, or for instance [24]), i.e., by 
definition: 



«(x) := -\ ^W^) f = -=*_ (3.54) 

2 A V<t(2/) Jo ^/ [l - y 2 ){\ - my 2 ) 



We have: 



w(a) = 0, u(b) = K(m), 
u(c) = K{m) + iK'(m), u(d) = iK'{m). (3.55) 



We have: 



ic-a sn(w, m) cn(u, m) dn(w, m) 

V a(x) = -i[d — a)(b — a) J - — - b _ a — . (3.56) 

\d — b (1 + ^sn 2 (u, m)Y 

Let us define such that: 

/ d-b 

Uqo := i / — , (3.57) 

Jo ^(l + y 2 )(l+my 2 y 



i.e. 



d — b Id — a (d — a) 

sn(w 00 ,m) = id- , , cn(M 00 ,m) = W- , , dn{u 00 ,m) = * -. (3.58) 

V o — a V o — a y(c — a) 

Then we define x : 

x = d + iy/(c - a)(d - b) [E^^m) - ( 1 - JTTT^ ] • (3.59) 



It satisfies: 



thus we have: 



/ X X ° dx = 0, (3.60) 

Jb a/ (x — a){x — b) (x — c)(x — d) 



Q(x) = X X ° (3.61) 

y/{x — a) (x — b) (x — c) (x — d) 



u^iuoc 9i((u(x) +u !X )/2K) 

A t = p KK' 

1 ] 6 1 {{u{x)-u O0 )/2K) 



(3.62) 



7 = ±y/(d-b)(c-a)e-*& ^ (3.63) 
4a e^Uoo/K, t) 
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F( n e 1 (u(x) + u(Q)e 1 (2 Uoo ) 

= 2 x ° (3 ' 65) 

we have the asymptotics fH\ 14]: 
d(T c %) - a(T c f ) - c(T c %) + 6(T C ^) <T c f ) - a(T c f ) + C (T C %) - 6(T C %) 

A — — /I 



(3.66) 

d(T c %) + a(T c §) - c(T c f ) + 6(T c f ) d(T c §) + a(T c f ) + c(T c f ) - fc(T c §) 



(3.67) 

Therefore, we shall now study W(x,T, V) in different regimes. 

4 The Birth of a cut critical point 

Let us choose the potential V and the temperature T c such that: 

• for T < T c we are in a one- cut case, 

W(x, T) = - (V'(x) - M_ (x, T) ^(x-a){x - &)) (4.68) 
with a(T) < 6(T) and 

M_(x,T c ) = (x-e) 2 " _1 Q(x), (4.69) 

where z/ > 1 is an integer, and Q(x) is a real polynomial whose properties are 
described bellow. 

• for T > T c we are in a two-cuts case, 

W(x,T) = - (v'(x) -M + (x,T)y/(x-a)(x-b)(x-c)(x-ctf) (4.70) 
with a(T) < b(T) < c(T) < d{T) and 

c(T c ) = d(T c ) = e , M + (x,T c ) = (x - e) 2u ~ 2 Q(x) , (4.71) 

• at T = T c one cut has vanishing size c(T c ) = d(T c ). With no loss of generality, 
we can assume that: 

a(T c ) = -2 , b(T c )=2, (4.72) 

and we write: 

e(T c ) = 2 cosh0 e = c(T c ) = d(T c ) (4.73) 
11 



The polynomial Q(x) must have the following properties: 

• degQ = d — 2v with d odd and d > 2u, 

• The leading coefficient of Q is positive, 

• Q has an odd number of zeroes in ]2, e[, 

• Q{x) < in [-2,2], 

• Q(e) > 0, 



\/x < -2, / Q(x)(x - e) 2 "-Vx 2 -4dx > (4.74) 

J X 

Vx>2, x^e, ^ Q(x)(x-e) 2 "~Va; 2 -4dx >0 (4.75) 

Q(x)(x - e) 2 "-Vi 2 - 4rfx = (4.76) 

V'(x) = Pol ( (x - e) 2 ^ 1 g(x)v / ^ 3 4 N ) (4.77) 

x— >oo \ / 



1 



T c = - Res (x - e) 2v - 1 Q(x)y/x 2 -Adx (4.78) 

2 oo 

Remark: notice that for all v > 1, it is possible to find a potential V(a;) and a tem- 
perature T c with such properties. Indeed, choose e and Q(x) with the above properties 
and determine V'(x) and T c bv 14.771 and 14.781 Notice also that it is always possible to 
find a polynomial Q(x) which satisfies the above mentioned conditions, indeed consider 
any real e > 2, and any real polynomial Q(x), of even degree d — 2v — 1, with positive 
leading coefficient, and with no real zero, then set: 

^xQ(x)(x-e) 2u - 1 Vx^Adx 
e = — — 7 (A 79) 

J 2 e Q(x) (x - e^Vx^Adx 

clearly, e e]2, e[, and then set: 

Q(x) ={x- e)Q(x) (4.80) 
In particular, one may choose d = 2u + 1 and Q = 1. 
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4.1 Example v = 1 

Let e > 2 be fixed. We write e = 2 cosh e . 
We consider the following quartic potential: 



V'(x) = (x 3 -(e + e)x 2 + (ee-2)x + 2(e + e)) , T c = 1 + ee (4.81) 



where e is given by f^(x — e)(x — e)V^ 2 — 4 = 0, i.e. : 

_ _ 0e cosh e - | sinh e (2 + cosh 2 e ) 
isinh0 e cosh0 e (5 — 2 cosh 4> e ) — 4> e 

4.2 At the critical point T = T c : 

At T = T c , both formula T4.68I and 14.701 reduce to: 

W{x, T c ) = ^ (v'(x) -{x- ef^Q^Vx^t) (4.83) 
which would correspond to an average large N eigenvalue density in [—2, 2]: 

p{x) = -^—(x - e) 2v ' 1 Q{x)V^x~ 1 (4.84) 
2irl c 



and one would have: 

7iv ~ 1 , Pn ~ (4.85) 

However, this is wrong, because the semiclassical asymptotics 13.221 are valid only if 
T 7^ T c , they break down at T = T c . It is the purpose of section to determine the 
asymptotic behavior of 7 n and (3 n near n = N and T = T c . For the moment, let us 
consider the limits of 14.681 and 14.701 near T c . 

4.3 Variations near the critical point, T < T c (one-cut) 

Let us consider the limit of 14. 681 near T c . Write T = T c + t, and t < 0, and: 

W(x, T) = - (V'(x) - M_ (x, T) v/(x-a)(x-6)) (4.86) 
At T = T c we have 

a(T c ) = -2 , b(T c ) = 2 , M_(x,T c ) = (a; - e) 2 ^ 1 Q(x) , (4.87) 

Then, make use of formula 13.291 and 13.481 i.e. 

_ 1 dM_(x,T) I % I j| 1 

2 <9T + 4 1 ' ' {x - a) + 4 1 ' J (x - 6) (x - o)(x - 6) 1 J 
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matching the pole at x = a gives: 

da 4 



dT (o-6)M_(o,T) {a - e) 2 "- 1 Q{a) 
which is finite at T = T c , thus, we find that to first order in t: 

t . _ t 



(4.89) 



(4.90) 



(2 + e) 2 "- 1 Q(-2) ' (e-2)2"-iQ(2) 
(notice that Q(-2) < and Q(2) < 0). Relation EOK1 implies: 

7~l+0(t) (4.91) 



Then, 14.881 reduces to: 



1 dM_(x,T) M_(x,T) - M_(a,T) M_(x,T) - M_(b,T) 

2 <9T M_(a,T)(a-b)(x-a) M_(b,T)(b - a)(x - b) { ' ' 

which is finite at T = T c , thus one gets the asymptotics of M_: 

<9M_(x,T) (x-e) 2v -\Q(x) - Q{a)) + {{x - e) 2 ^ 1 - (a - e) 2 ^- 1 )Q(a) 



2 



dT (a - e) 2 "- 1 Q(a) (x - a) 

(x - ef v -\Q{x) - Q(b)) + ((x - e) 2 "" 1 - (6 - e) 2 "- 1 )g(6) 



(6 - e) 2 *'" 1 Q{b) (x-b) 
{x - e) 2 "' 1 Q{x) - Q(a) (x - e) 2u ~ l - (a - e) 2v ~ l 



(a - e) 2 "- 1 Q(a) x-a (x - a) (a - e) 21 " 1 

(x-e) 2 "- 1 Q(x)-Q(b) (x - e) 2 "" 1 - (6 - e) 2 "" 1 



+ 



(6 - e) 2 "" 1 Q(6) x-6 (sc - b) (b - e) 2 ^ 1 



(4.93) 

in particular in the vicinity of x = e one has: 



M_(x,T) ~ (x-e)^ _1 Q(x 

r2i/-2 

+ 2 



- e) fe ((2 - e)-* -1 - (-2 - e) - * -1 ) + 0(x - e) 2u ~ 

k=0 

(4.94) 

Note that the zeroes of M_(x,T) in the vicinity of e, are the 2v — 1 th roots of unity: 

/ Of \ V^-i 
M_(x,T c + t) = <-> x = e+ ^____j +0(t 2 ^- 1 ) (4.95) 

Using I33TI we get: 

1 t / 1 1 

7n ~ 1 - 7 T- ono„_i^/on + 



4 V (e - 2) 2 "- 1 Q(2) (e + 2) 2 "" 1 Q(-2) 



^ 2 V (e-2) 2 "- 1 Q(2) (e + 2) 2 --!g(-2) 
where n = N(l + t/T c ) (4.96) 
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4.4 Variations near the critical point, T > T c (two cuts) 

Let us consider the limit of l4.7Ul near T c . Write T = T c + 1, and t > 0,and 

W(x, T) = - (V\x) - M+(x, T) y/(x-a){x-b){x-c){x-dj} (4.97) 
At T = T c we have 

a(T c ) = -2 , 6(T C ) = 2 
c(T c ) = d(T c ) = e = 2cosh0 e , M+(x, T c ) = (x - e) 2 ^ 2 Q(x) , (4.98) 

and at T > T c , a + 2, b — 2, c — e, d — e, and M + (x) — (x — e) 2u ~ 2 Q(x) are small. In 
particular, we write: 

M+(x, T) = H(x, T)Q(x, T) (4.99) 

where H(x,T) is a monic polynomial of degree 2v — 2 which contains all the roots of 
M + close to e, and Q(x, T) is the remaining part. In other words, H(x, T) — (x — e) 2u ~ 2 
is small and Q(x, T) — Q(x) is small in the small t limit. 

We use the notations of section 13*31 The biratio is thus: 

(b — a)(d — c) 4 ., . d — c 

m = ) r(d-c) s— 4.100 

(c-o)(d-6) e 2 -4 v ; sinh 2 e V ' 

we see that we have to consider the limit m — > 0. In that limit 13.571 becomes 



' d-b / e-2 

b - a ay . / V 4 fly 



u 00 = i = ~j / * =z^ (4.101) 



E(uoo) -Uoo = im 



1 d=b „ 
!>-<• y ay 



o v/(l+y 2 )(l + my 2 ) 



^i=2 r, 
4 



(4.102) 



y dy sinh < 
im / — , = im 







And, as can be found in any handbook of classical functions [], we have the small m 
behavior: 

m ~ ~ (4-103) 
K'[m) mm 

Since for small m one has \r^— I 3> m. 13.591 becomes: 

I mm I ' 1 1 

20 e sinh0 e 

ox := x — d ~ (4.104) 

mm 
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Notice that: 



x — c = xq — d+d — c^xq — d + m sinh 2 <p e ~ x — d ~ 5xq (4.105) 



From l3~2*9~l and 13". 6 1 1 we have: 

1 dM+(x,T) 1 M+(x,T) da{x,T) x - x 



2 dT 4 <r(x,T) dT a{x,T) 

Matching the pole at x = c gives: 

dc —A(d — c + 5x ) — A5xq 1 



(4.106) 



dt (c- a)(c-b)(c- d)M+(c,T c + t) (e 2 - A)Q(e) (c - d) H(c) 

and matching the pole at x = d gives: 

dd —45xo A5xq 1 

di = (d-a)(d-b)(d-c)M+(d) ~ ~(e 2 - 4)Q(e) (d - c) #(d) 

and matching the poles close to e gives: 

25x 1 (H(x)-H{d) H(x)-H(c) 



(4.107) 



(4.108) 



^ (e 2 -4)Q(e)d-cUi-d)5(d) (ac - c)H(c 



(4.109) 



The following guess solves the 3 equations 14.1071 14.1081 14.1091 to small t leading 
order: 

j_ j_ 



where ( is a positive real number, and G is a degree 2z/ — 2 even monic polynomial, 
which will be determined below. 

For later convenience, we also define the following positive constant: 

sinh0 e Q(e) ^ 
Using ansatz l4*. 1101 and I4.104[ we have in that limit: 

4W> e sinh0 e f A -iio\ 

ox — (4.113) 

hit 

Then, inserting- KTW\ and l4~TTTl into KW7\ and 14.1081 we get: 

AC 2 = — — = — — (4.114) 
^ G(-2C) G(2() { ] 

16 



Then, setting x = e + £ (— j-M 2 " , and inserting 14 . 1 1 1 1 into 14 . 1 091 we get the following 
equation for G: 

(2u 2) G(() CC'(C)- C ( G ®~ G W G(0-G(-20\ 

(2, - 2) G(0 - (0 - ^ ^ (e _ 2C)G!(2C) - (e + 2C)G( _ 2C )J ( 4 - 115 ) 

which using |4*~114I becomes: 



> -2) 0(0-^(0 = ^53^5 (G(0"G(20) (4.116) 



the solution of which is: 



G(Q = V ^- ( 2k ^-i-x) = Pol J (4.117) 



or: 



G(2Ccosh^)= C 2l - 2 y: f 2 "" 1 ! Sinh(2j u +1) ^ (4.118 
G(2C) 1 2z/! C 



3=0 

In particular 



^-2 2i/-l!i/l 4C 2y 
i.e. the parameter £ is determined by: 



(4.119) 



^ V 2 2z/! / Vsinh0 e g(e) 2i/! 7 1 ' ' 



In that scaling regime, we have: 



i 



m (-7^- J (4.121) 

sinh 2 (/) P V lnt ' 



i.e. this corresponds to a torus of modulus 



and, using EiESJ 

7 



X' —i —i 
r = i mm lnt (4.122) 



—-\/ (a — b)(c — a) e 



4^ vv yv ' 6 1 (u 00 /K,r) 

1 # 



lnm 

~ 1 

lnt 

(4.123) 
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we also find that the filling fraction in the [c, d] cut is of order: 



1 <- d 



2nT CJ , 



p(x) dx 
t sinh0 e Q(e) f K 



In t 7rT r 



J-2C 



-2C 

t smh0 e Q( e ) r^^GioV^e^e 



(4.124) 

then, integrating by parts and using 14.1161 we find: 



T c lnt 


2vk 


t 


sinh0 e Q(e) 


V c lnt 


2vk 


t 


sinh0 e Q(e) 



L-^^fe G(2C)df 



•2 



* de 



2C v 7 ? 3 ^ 



(4.125) 



T r lnt 2i/ 



C 



i.e. 



e (T c + t)~-—L_2u<t) e (4.126) 
J c lnt 

This means that for n = N(l + t/T c ), the average number of eigenvalues located near 
e is: 

n — N 

k ^~kjr 2u<pe (4127) 

i.e. the eigenvalues start to explore the potential well near e when n — N ~ IniV. 
According to 13.661 and 13.671 the coefficients 7„ and /5 n vary between: 

i 

1 + C 2 " < 7n < cosh0 e (4.128) 

2C (ln^)" ^^ e " 2 ( 4129 ) 
where n = N(l+t/T c ) 
The transition takes place on a scale of order p ~ In A/". 

4.5 Order of the transition 

From l3~3"4l and l4*.91l we have below T c , i.e. for t < 0: 

()iF (T c + t) = -2 In 7 ~ t (——*— + ___L__ > ) (4.130) 



<9t 2 v L y ' 2 V (e-2) 2 "- 1 Q(2) (e + 2) 2 ^ 1 Q(-2) 

1 



and above T c , i.e. for t > 0, we have from 14.12*31 

^CT. + «) = -^ 7 ~^ (4.13!) 



The second derivative of the free energy is continuous, but the third derivative is not. 
therefore we have a third order transition, with logarithmic divergency, cf figEl 



2 ; 




9F 








3t 






t 







Figure 3: Behavior of the second derivative of the classical free energy with respect to 


t near the transition. 



5 Mean field asymptotics for the partition function 

We compute the partition function for a potential: 

V ro (x) = V(x) +r ln(x - x) (5.132) 
where r is assumed of order 1/N: 

T c 

r ° = _a i7 (5.133) 

5.1 Mean Field theory 

We use the same idea as in [3] . We split the 2-cuts integral into 1-cut integrals. Let us 
say that there are k eigenvalues in the new cut (near e), and n — k in the old cut [a, b\: 

Tl 

n\ Z n (—T c ,V ro ) 
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n / \ „ k 

~ El*) / da; 1 ...dx fc A 2 (x i ) n^o-^re"^^ 
fe=0 \ ' »=i 

(n - fc)! Z n _ fc (r c , V ro ,rJ 

iv 

(5.134) 

where 

3=0 



(5.135) 



and Zn-j-iTc 1 ^, V rj )is a one-cut integral: 



^- fe (T c ^— -,V rj ) := -J— / dx^.-.d^A 2 ^) TT e-€^W 



= # n „ fc e ^-^^'^ (5.i3 6 ) 
That gives: 

Zn(jjT c , Vr ) 



~ do; 1 ...dx fc A 2 (x 4 )n(xo-x i re-€ y ^e ^' 

(5.137) 



In other words, we integrate out n — k eigenvalues, and consider the integral over k 
eigenvalues only. The k remaining eigenvalues are submitted to the potential V, as well 
as their mutual Coulomb repulsion, and the mean field of the exterior n — k eigenvalues. 

Since F n (T,V rj ) corresponds to a one-cut distribution, it can be evaluated with 
standard semiclassical technique (see sectional), and in particular, it has a large n 
expansion: 

F n (T, V rj ) ~ F(T, V rj ) + -F (1/2) (T, V r ,) + — F (1) (T, V r .) + O(-) (5.138) 

and each term of the expansion is analytical in V r . . 

We want to evaluate F(T, V r ) in a regime where T — T c is "small" (we make that 
more precise below) and tq and the r/s are of order 0(1/N). 

we first do a Taylor expansion in T — T c and the r^'s: 

k ffp 1 f) 2 ~~F k 2 



dr,- 2^ drjdrj N 2 ' 

i=0 ij J 

d 2 F 
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r)F 1 8 2 F 1 

+ (T -T C )— + -(T- T c f— + o(^) (5.139) 

where all derivatives are computed at T = T c and = (see section l3~2l and EPj) . 
That expansion is valid only if k <C iV. 
Thus we have: 

-1- Z n ( ™ T c , K ) ~ e"* F(Tc ' V^,vV^o e -^o.o 
if n A 7 " 



it 



1 H n ^ k c _^j n _ N _ k)FTca(n _ N _ k) F Trnc _ ^-N-kf TrT 



k\ K, 

k=0 ' 



/ d Xl ...dx k A 2 ( Xi )Y[\ 

Jxi>e 



_ N 

Xo-Xife 



(5.140) 



k _ _ _ _ 

J]V£^; e 2(n-N-k)F T , rj e -2aF ro , rj j T e ~2F ri , rj 



Notice that: 



~ (2vr)- fc (l - -)-Vi2 (1 + 0(1/n)) (5 . 141) 
5.2 Computation of the derivatives 

Now, use formula given in section I3~21 and l3~3~l in the one-cut case, and get (derivatives 
taken at T = T c and T{ = 0, and take into account that V e s(b) = Kff(e)): 

—F = --(V eS (x i )-V(x i )) , —F = V eS {e) (5.142) 
^ -F = In — X \ - % . , ^F = -lnA'(x i ) (5.143) 



dridrj — ^( x j) dr. 

d 2 ^ , . , , <9 2 ^ 



F = lnA(x i ) , ^^ = (5.144) 



<9T<9r, v " ' <9T 2 

Moreover: 

9T{T, V) a + b dT(T, V) 1 

-^^ = — = ' ^^ = A^y (5 - 145) 

The effective potential behaves in the vicinity of e as: 

V eS (x) ~ V cS (e) + V -^f>- (x - e) 2v (5.146) 

x^e lv\ 



V£ v \e) = 2sinh0 e Q(e) 
2z/! 2f 



(5.147) 
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In the limit where the CC{ S £1X6 close to e, we have: 

1 



and 



A(£j) — A(Xj) Xi,Xj^e A'(Xj) 

A(xi) ~ e^ e 



2 sinh e e 



(5.148) 



(5.149) 



5.3 Result 

write n = N + p: 

Zn+ p ( n P T c , V ro ) ~ 



f e* e 
\ 2 sinh f 



e -^(Vcft(^0)-V(x o )) 



H N+p _ k e _N_ (p _ k)Vcs(e) e a(p-fc)^ e 

^ fc! 



k=0 



[ d Xl ... dx k A 2 (x t ) f[(x - x t ) a e~% v ^ 

k , 
JJ e -^(^«ff(*i)-n*i)) e 2fc(p-fc)* e j 



i=l 



\ 2 sinh e 



2afc+2fc 2 



e -% F ( T - y ) e -F (1) (T c ,y) g-^CVrffCxoJ-VCxo)) e -pf V cH (e) 



N +P (0~\p-k f 
' (2fc+a)(p-fc)</> e / _ 



fc=0 



k\ 



\ 2 sinh <f) ( 



2(fc+§) 2 



/ dxi . . . da: fc A 2 (^) JJ(a; -x i ) a e" 



;(V r eff(x i )-y e fr(e)) 



(5.150) 



Then we rescale: 

Xi = e + 
and we get: 



/2sinh0 e Q(e)\ 2v -i / 2 sinh</> e Q(e)\ 2 " 

' ' Vl , x = e + N- I f^ 1 ] y (5.151) 



,N + V . aN 

Z N+p (—^T c ,V ro )e-^ 



V(x ) 



N 



Hn e -^ F ^ e -F W (T c ,V) e -( P+ a/2)£v cff (e) Q -a^ 

(2sinh0 e )^ (2vr) p 

N+p 2 

jV~ (fc+2 ^ /2) e ( fc + a /2)2p0 e e (fc+a/2)a<^ e ^-(fc 2 +afc) ( 27r )- fe 

fc=0 
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^ J dj/i . . . dy fc A 2 (^) - Vl ) a e 



k 

I chh r\ih. A 2 (nA 

i=l 

(l + 0(i\T&)) (5.152) 



where we have defined: 



J 4:=(2 3 mh^( 23in "^ Q(e) ) (5.153) 



Eq. (|5.152j1 is valid only up to 0(N~^) because Eq. ()5.146j) . Eq. (j5.148j) and 
Eq. (|5.149jl are valid only to that order in the regime of Eq. (j5.151|) . 

5.4 The effective matrix model 

Let us define the matrix model in the potential 

Let (k,u be the partition function of the k x k matrix model in the potential ^— : 

1 f x 2v 
(k,v--=-jj / d Xl ...dx k A 2 ( Xi ) Y[e~^ (5.154) 

Notice that for v — 1, this is the gaussian matrix model, and we have: 

In C*,i = |ln27r + In (JJj!) = \ n H k + ^\nk - ^k 2 (5.155) 

3=0 

We define the amplitude: 

A k := A- k2 (2n)- k ( k ,u (5.156) 

We also introduce: 

hk v := = 2irA 2k+1 ^±1 (5.157) 

Qk,u A k 

^ := V^ = X V^T" (5 ' 158) 

and the associated orthogonal polynomials 

P k (y) := /^•••^fo) 2 nUv-Vi)^ (5 15Q) 

fdyt... dy k A(yi) 2 [jJU e"^ 



:= d-r^-PMe- 3 ^ (5.160) 



and their Hilbert transforms: 



J dx 1 ... dx k A(xi) 2 n?=i e 2 " 



:= " , > (5.161) 

/ dxi . . . dx k A(xi) 2 Yl i=1 e ^ 

4>hM '■= \[hk~vPk{y)z^ (5.162) 
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5.5 Partition function 



Thus we have for a = 0: 



Zn+ p { j-f P T C , V) 



,_ N+p 
-F(T C ,V) f W (T c ,V) 



H N e -^ F{T - v \-F' \t c ,v) e - P f v cff (e) (2?r)P Arfe e 2k ^A k {\ + O(AT^)) 



(5.163) 



k=0 



(eS p+1 ^ e2We e2fc0e 4*) (est 1 e2We e_2fc</,e ^) 

(ES rie * ( A k) 

(1 + 0(JV - ^)) (5.165) 

5.6 Orthogonal polynomial 

According to Heine's formula (cf Eq. (|2.12j) ). we have: 

, Z n (T e ±V(x)-%\n(Z-x)) 
Z n (T c jf, V) Z n+ i(T c ^, V) 



n _l(O = . e 2T V ^> 



ip N+p {x Q ) 



Z n {T c %,V)Z n _ x {T c ^,V) 

Zn+ p {^^-T C) V ro ) e 

Zn+ p {^tT c , V) Z N+p+ i( N+ £ +1 T c , V) 
T 

r = ~ , a = 1 (5.166) 



thus, in the regime: 
using Eq. (|5.152|) with a = 1 we get: 



N 

i 4 sinh 2 e 

x = e + iV~^ -—^y (5.167) 



A 



2 sinh e 

YSSn-**^ eg+VW; e( fc+1 / 2 )^ V3j^T VyQ/) 



(ES P+1 ^ e2fcp0e e2fc</,e 40 (£2f ^ e2fcp0e 4,) 
(1 + 0(JV - ^)) (5.168) 
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5.7 Hilbert transforms 

According to Eq. ([2. 15)1 . we have: 

(f) N+p -i(x ) = 



Z N+p (^T c ,V ro )e^) 



N 

thus, in the regime: 

i 4 sinh 2 



Zn+ p -i( N+ n T c , V) Z N+p (^^T c , V) 
T 

r = , a = -1 (5.169) 



x = e + N-^ T-^V (5.170) 



using Eq. (|5.152|) with a = — 1 we get: 



4> n+p -i(xq) ~ N*» 



.4 



2 sinh e 

V (ES ^ e 2 *^ ^ e2fcp * e e " 2fc</,e ^) 

(1 + 0(AT&)) (5.171) 

Notice that shifting k — > + 1 we have: 



07V +p _l(x o ) ~ iVs- 



.4 



2 sinh e 



^2 



'(E2? ^"^ e 2 *^ A A )(X;S ,_1 JV"£ e 2 ^ e" 2 ** A fc ) 

(5.172) 



5.8 Computation of (3n+ p 

We start from Eq. (j2~TTj) 



iV_,n I / , ^ %) . 2 



(£x, + ^T„_ fc (T c ^, V)) e-S'-f-^ (5.173) 



We have: 



T n _ fc (T,V) ~ T(T c ,r) + (T-T c ) — (T^^ + ^r, — (T c ,^) + 0(l/iV 2 ) 

i 3 



25 



IT 

' c ^A(x,)- 1 + 0(l/iV 2 ) (5.174) 



N 

3 



Therefore: 



— T(-T V) ~ - / dxi dx t e-€^)A 2 f 

T ^ jV ZfiT fc' / - ' ' ' 



Xi 



N .F n _ k (T c Hj±,V) 



^ c ' y fc=0 

fc 

(E-o-^, ; ) <• 

3=1 

2sinh ^ V fc / dx, dx, e~^ v ^A 2 (x-) 

_^F„_ fc (T c ^,V) 

e c 

2sinh0 e ^^ 7V 2 (5.175) 



and: 



6 Asymptotic regimes 

Consider p of order In AT: 

it 

p= -In AT , u finite. (6.177) 

6.1 Possible asymptotic regimes for the partition function 

Then Eq. f!5.163|) becomes: 



2 — N +P 

fc=0 



(6.178) 

It is clear that the sum over k is dominated by the values of k for which the exponent 
of A" is maximal, i.e. for which 2mA; — k 2 is maximal. This means that for u < 0, the 
sum is dominated by the vicinity of k = 0, and for u > 0, the sum is dominated by the 
vicinity of k = u. The sum is then well approximated by a few largest terms. 
Let us denote u the positive integer closest to u: 

«:=[« + 1/2] if.ao 
it := it it < 
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where [.] denotes the integer part. 
Define also: 

e u := sgn(w - u) if u > , . 

e u := 1 if u < (6 ' 180) 

We always have: 

e u ( M -M)<i (6.181) 

The largest value of 2uk — k 2 is obtained for k — u, and the second largest value is 
obtained for k = u + e u . The difference is: 

(2uu - u 2 ) - (2u(u + e u ) - (u + e u ) 2 ) = 1 - 2e u (u - u) (6.182) 



12 3 4 




Figure 4: Behavior of e u (u — u 



2' 



Remember that our asymptotics for Z n are valid only up to order 0(N 1//2v , i.e. 
we want to have: 

e u {u-u)-\>-\ (6.183) 
which implies that our asymptotics for Z n are valid only if u > and u N. 

6.2 Possible asymptotic regimes for the orthogonal polynomi- 
als 

We have similar considerations for the asymptotics of orthogonal polynomials, excpet 
that the sum over k is now shifted by 1/2. This gives different regimes. 
Eq. ()5.168j) becomes in that regime 



ip N+p (x ) ~ iVs. 



A 
2 sinh ( 
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(Ek=o +1 N ~£ e 2fcp * e ^ A k ) (Zk=o e 2kp * e A k) 

(1 + 0(JV~^)) (6.184) 

The sum in the numerator is dominated by the largest values of 

2(k + l/2)u - (k + 1/2) 2 , (6.185) 

i.e. by k = [u] = u + 

The second largest term is obtained for k = u — Notice that if u > |, the two 
largest terms are always u and u — 1. The difference 

(2{u + l/2)u - (u + 1/2) 2 ) - (2{u - l/2)u - (u - 1/2) 2 ) = 2{u - u) (6.186) 

Since our asymptotics are valid only up to order O(N^), the subleading term 
should be discarded if \u — u\ > |, i.e. if it < ~ or if w is half-integer. 

This implies that our asymptotics for ip n are valid only if u > | and w ^ N + ~. 

6.3 Asymptotics in the regime it > 0, and u not integer or 
half- integer 

From now on, we write: 

u 

n = N + p , p = IniV. (6.187) 

2z/0 e 

In this section, we assume that u > and u not integer or half integer. The sum 
over k in Eq. (|6.178|) is dominated by the terms k = u and k = u + e u . 

6.3.1 coefficient j n 

Thus we have: 



(6.188) 

We also obtain: 



Ok-TT—iT 2 / W — U\ — 7T 1 



ki — U I — -i /I — . 



//.v-,, ~ 2 7T e-^ Kff{e) e 2 ^ ^1 + 2e u e e ^ e sinh e iV^^ + 0(N~^) 

(6.189) 
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and: 



Jn+p ~ 1 + 2 sinh <f> e N 
(6.190) 

In AT 



2. s±u+e v 



+ 0(N~~> 



jN+p is nearly periodic, with period ^-r-- The amplitude is minimal of order N ^ for 
u integer, and is maximal of order 1 for u half-integer, cf fig. 16.3.11 



Yn 




InN ,--■"'"7*' 


ch(p e - 






-l/2v 
N 


1 - 












N 


n 




Figure 5: Behavior of 7„. 



6.3.2 coefficient (3 n 

The sum in Eq. (I5.176J1 is dominated by k = u and k = u + e„, i.e. 

(3 N+P ~ 4 sinh 2 e jV 1 ^^ 1 e £ ^ e ( " +£ " } 



AjT 



6.3.3 Orthogonal polynomials 



(6.191) 



The two largest terms in the numerator of Eq. ()5.168|) correspond to k = u and k = it— 1 
(not necessarily in this order), thus: 



1pN+ P (x ) 



A 



2 sinh ( 



N U ~^ ^ u{y) + e -^/2 ^ ^_ lAy) 



1 + cosh 6 e V" 2 e e ^ e — 



(i + o(j\r*0) 



(6.192) 



we also have: 

^Ar +p _i(a;o) 



.4 



2 sinh ( 
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e"^ 2 ^ Ay) + e ^/2 V%_ v (</) 



1 + cosh (p e N~ 



(1 + 0(JV - ^)) 



(6.193) 



6.3.4 Hilbert transforms 



Similarly, using Eq. (J5.172)) . we find that the Hilbert transform of the orthogonal 
polynomial 7r n , are asymptoticaly given by: 



A 



2 sinh (j) e 



e"^ 2 fajy) + e^ 2 J+p ^ v {y) 



2\u — u\ — l 



1 + cosh <p e N •*<> e 



(6.194) 



and 



4>N+p{Xo) 



A 



2 sinh e 



7V^? e^ 2 <hAv) + N- U ^ e-^/ 2 fa^y) 



1 + cosh e N 



2\u — u\— 1 



2v e u 



(6.195) 



ii 



6.3.5 Matrix form 

The matrix: 



V'n-l(^) 0n-l(^) 



is in that regime (w > 0): 



2sinh0 e \e~1^ e e^ 



J V ^My) <My) 



(6.196) 



(6.197) 



where 



R — di&KlN 2» \l — : — ,Ni" 



L = 1 + cosh e iV~ 



. 2|u— TT| — 1 t4" 



,1,. I 



(6.198) 
(6.199) 
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6.3.6 Kernel 



The kernel K n (x,x') is given by the Christ ioffel-Darboux formula: 

x — X 

We find: 

K n {x, x >) ~ N ^ 2 A 7 „ Mv)^M-**-iAy)My') (6201) 

4 sinlr e ' (y - y') 

i.e. 

K n (x,x')~ K %v (y,y') £ (6.202) 



6.3.7 Large it limit 



Using Stirling's formula (cf appendix H), we find that for large u, all those asymptotics 
match with the classical limit of section 4.4. 



7 Conclusion 

We have computed the asymptotics of orthogonal polynomials in the birth of a cut 
critical limit. This corresponds to the appearance of a new connected component for 
the support of the eigenvalue density, away from other cuts. 

We have found some universal behaviour, which depends only on the degree v of 
vanishing of the density at the new cut, and 2cosh0 e which parametrizes the distance 
between the new cut and the old cut. The parametrix near the new cut, is simply 
the system corresponding to a model matrix model in the potential x 2u , and with 
u — [| + 2z/0 e ^^] eigenvalues. 

This new universal behaviour does not seem to correspond to a conformal field 
theory (unlike previously known critical behaviours), because the exponent of iV is not 
constant. 

It would be interesting to complete the "physicist's proof" presented here, with a 
mathematical one, for instance using Riemann-Hilbert metods as in |2| lllj. 
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Appendix H Stirling formula and other asymp- 
totics 



Stirling formula: 

1 



n\ ~ n n e- n V2^(l + — + . . .) (H.l) 



from which we deduce: 



1 3 n 1 

ln(l...(n- 1)!) ~ -n 2 \nn- -n 2 + -ln27r - — Inn + 0(1) (H.2) 

and 

Inn . 

\nH n ~ nln27r - — + . . . (H.3) 

Asymptotics of (k,v For large k we have: 

lnCfc*, ~ — lnife- — k 2 + -lnk + 0(k) (H.4) 

Appendix I Elliptical functions 

We introduce a few definitions about elliptical functions |24 : The elliptical sine func- 
tion sn(w, m) is defined by the following identity: 

, snM d 

u= / . (LI) 

Jo y/(l-y 2 )(l-my*) 

The complete integrals are defined by: 

K{m) := C dy = (1.2) 

Jo V(l-y 2 )(l-my 2 ) 

g'( m ) ■= y(i - m ) = I" , dy (1-3) 



Wu,m) l_ mv 2 
E( Uj m):= j( J-—JL d y (1.4) 



E ( m )-= I \l 1 ™y2 d y > ^'H := E{1 - m) (1.5) 



When m — > one has: 

7r / m 9m 2 _ . o N 

if ~ -IH 1 h • • • + 0(m s 

2 V 4 64 v ' 

... , 1 /, m 9m 2 

if 7 ~ In — = 1 + — + — — + . . . + O m 3 
y/m \ 4 64 
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E 



E 




+ 0{m 2 ) 



+ ... 



+ 0(m 3 ) 



(1.6) 
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